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ABSTRACT 

The ultra-luminous X-ray source (ULX) M82 X-l is one of the most promising intermediate mass 
black hole candidates in the local universe based on its high X-ray luminosities (10 40 — 10 41 erg s _1 ) and 
quasi-periodic oscillations, and is possibly associated with a radio flare source. In this work, applying 
the sub-pixel technique to the 120 ks Chandra observation (ID: 10543) of M82 X-l, we split M82 X-l 
into two sources separated by 1.1". The secondary source is not detected in other M82 observations. 
The radio flare source is found to associate not with M82 X-l, but with the nearby transient source SI 
with an outburst luminosity of ~ 10 39 erg s _1 . With X-ray outburst and radio flare activities analogous 
to the recently discovered micro-quasar in M31, SI is likely to be a micro-quasar hidden in the shadow 
of M82 X-l. 

Subject headings: Galaxy: individual(M 82) — X-rays: binaries 


1. INTRODUCTION 


with X-ray luminosities of 


10 4 ° erg s 1 and quasi-periodical oscillation (QPO) 


The ULX M82 X-l, 

> 

behavior, is one of the most promising interme¬ 
diate mass black hole (IMB H ) can di dates in the 
local universe (ICollura et al.l 1994 Griffiths et al.l 


20001: Matsumoto et al.1 20011 : IStrohmaver et all 120031 : 


Pasham et al. 2014f) . The discovery of X-ray spectral 


state transition in M82 X-l shows it has X-ray luminosi¬ 
ties of Lx ~ 10 41 erg s _1 at the thermal-dominant 

state, suggesting an IMBH with mass ~ 200 to 8OOM 0 
(|Feng fc Kaaretl 2010t). The d iscovery of the 60-120 milli- 
Hz QPOs (IStrohmaver et all 120031: iFiorito fe Titarchukl 


2004 : Mucciarelli et al|200~ jCaballero-Garcfa et al .11201.4 


Pasham fc Strohmaved I2013T) suggests the strong X-ray 
emission originates from the accretion disk instead of a 
jet, and suggests an IMBH of ~ 2 x 1O 4 M 0 if these low- 
frequency QPOs corresponds to the Keplerian frequency 
of the innermost stable orbit. Pasham et al. (2014) re¬ 
analyzed 6-year RXTE X-ray observations and revealed 
the high-frequency, 3:2 ratio-ed, twin-peak QPOs of 3.3Hz 
and 5Hz, which, in combination with the low-frequency 
QPO revealed by an XMM-Newton observation, suggests 
an IMBH of 415 ± 63M 0 under the relativistic precession 
model (Motta et al. 2013). 

Many efforts have been made to search for the ra¬ 
dio counterpart of M82 X-l, and two radio sources were 
discovered at the vicinity of this ULX. A stable radio 
source 0.8" away from M82 X-l is proposed to be a 
supernova remnant (SNR) and not related to M82 X-l 
(ICollura et al. 1994 : McDonald et aJ]l2002t [Kording et al.l 


120051 : iGendre et ahl 201311 . The other radio source was de¬ 

tected (at R.A.=9:55:50.19, Dec=+69:40:46.0) through¬ 
out t he 7-hr 1981 VLA obser vation at a 9mJy flux 
level (|Kronberg fc Sramekl Il985h . but late r was not de¬ 
tected in the 2003 VLA observat ions Kording et all 
(|2005ll or the 2005 VLA observations dKaaret et al.l 2006( 1. 


iKording et al.l ( 2005!) re-examined the 1981 VLA data and 
found a 0.5" offset between the flare source and M82 X- 
1. Given the Chandra positional uncertainty of 0.6", it 
is not clear whether it is really associated with M82 X-l 
and suggesting relativistic beaming effects as causing the 
high X-ray luminosities. A careful work on astrometry is 
required to achieve a better understanding of M82 X-l. 

In this paper, we re-analyze the Chandra ACIS obser¬ 
vations with the sub-pixel technique to obtain improved 
spatial resolution and centroid accuracy. We describe data 
analysis and results in § 2 and discuss the implications of 
our work in § 3. 


2. DATA ANALYSIS & RESULTS 
2.1. Spatial Analysis 

M82 X-l has been observed 22 times by ACIS onboard 
Chandra between 1999 and 2014, with off-axis angles from 
below V to above 4'. The Obs IDs are: 361, 378, 379, 380, 
1302, 2933, 5644, 6097, 6361, 8190, 10025, 10026, 10027, 
10542, 10543, 10544, 10545, 10925, 11104, 11800, 13796 
and 15616. In this work we focus on Obs ID 10543 (PI: 
Strickland), which had the longest exposure of 120 ks with 
M82 X-l placed only ~ V from the aiming point. The 
energy range used in this paper is 0.3-8 keV unless stated 
otherwise. A careful examination of the Chandra image 
shows that the image size of M82 X-l is larger than other 
nearby sources (e.g., J095551.2, see Fig[T]). motivating our 
further analysis. 

We reprocess the Level 1 event files using the chan- 
dra_repro command in Chandra Interactive Analysis of 
Observations ( CIA 0,version 4.5) software. Following the 
standard pipeline^, we remove time intervals with sig¬ 
nificant flares, as well as bad columns/pixels and events 
of bad grades, to produce the Level 2 event files. The 
check-vf-pha parameter was set to no to avoid false- 
removing of good events. The sub_pixel parameter was 
set to EDSER to enable sub-pixel analysis. Fig |T| shows 
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the 1/8 pixel ACIS 0.3-8 keV image centered around M82 
X-l after reprocessing. Clearly, M82 X-l is split into 
two fuzzy but discrete sources, while other nearby ULXs 
are not. Applying wavdetect to the image results in an 
detection of two sources, as presented in Table [T] and 
FigdJ The brighter one is M82 X-l itself, and its position 
(RA=9:55:50.123, Dec=+69:40:46.54) is consistent with 
previously detected positions (e.g., Griffiths et al. 2000). 
The second source (hereafter SI), with RA=9:55:50.245 
and Dec=+69:40:45.67, is located 1.1" to the lower-left 
of M82 X-l. Not detected by wavdetect is extra photon 
concentration to the lower right of X-l, which is probably 
caused by the PSF feature like other nearby ULXs (e.g., 
see J095551.2 in FigOQ). In ObsID 10543, M82 X-l suffers 
from severe pile-up effects (see §2.2), which may affect its 
centroid. In comparison, the HRC observations do not suf¬ 
fer from pile-up effects and are expected to give accurate 
centroid positions for M82 X-l and other sources. We have 
analyzed the HRC observations and listed the centroid po¬ 
sitions for matched sources between the ACIS and HRC 
(Obs ID 8189) observations. As shown in Tabled similar 
displacement of ~ 0.6" exists between the ACIS and HRC 
observations for all matched pairs, including X-l despite 
the pile-up effects. 

As Chandra’s PSF could be strongly distorted and cause 
fake detection in +Z direction (Chandra instrument lay¬ 
out directions are used in this paper for convenience^, we 
need to ascertain whether they are real sources or just PSF 
artifacts. 

The relative position of SI strongly suggested SI is a real 
source. The relative position of SI to M82 X-l in CCD is 
— Y, which is different from the possible spurious source 
region (Fig [lj the cyan region centered at 15° counter 
clockwise to +Z) or the readout streak region (Fig [T]). 
Moreover, the distance from X-l to SI is ~1.1", greater 
than the PSF distortion region(~0.8"), making it unlikely 
to be a PSF artifact. Indeed, none of other nearby bright 
sources show photon concentration in — Y direction like 
SI, again supporting that SI is a real source. 

2.2. Spectral Analysis 

The spectra and the respective rmf and arf Hies of M82 
X-l and SI are extracted using specextract command in 
Cl AO with 0.7" circular source regions (Fig[T]). The back¬ 
ground region for X-l is a rectangular region, and that for 
SI is an annulus centered at SI, with inner radius 0.8" 
and outer radius 1.4". The X-l region is excluded when 
extracting spectrum for SI, and vice versa. The respec¬ 
tive count rates are roughly 0.029 c/s and 0.016 c/s in 0.3-8 
keV band for X-l and SI. The 0.7" extraction radius of 
SI may miss certain amount of hard photons outside, thus 
lead to a biased spectrum. Yet our results are qualita¬ 
tively valid due to the following reasons: i) the 0.7" radius 
is comparable to the local ~85% energy encircled radius 
as suggested by the PSF file, so the number of leaked hard 
photons is limited; ii) further correction on the spectrum 
will include more hard photons and thus a harder spec¬ 
trum. The spectra are then loaded to xspec and analyzed. 
We do not tend to focus on the spectra investigation, so 
an absorbed power law (pha*powerlaw) model is used as 
a characterization. 


As shown in Table [TJ the spectrum of M82 X-l has 
rm,x ~ —0.7, even h arder than the one reported by 
lAgrawal fe Misral (|2006l) . Such a hard spectrum is un¬ 
physical and likely caused by heavy pile-up effects, as in¬ 
dicated by the strong read-out streak. Because the read¬ 
out streak contains photons of M82 X-l before pile-up, its 
spectrum (with a 0.3-8 keV count rate of 0.009 c/s) is ana¬ 
lyzed as a proxy for M82 X-l. The spectrum of the readout 
streak can be well fitted with an absorbed power law with 
hstreak = 1-851 q33 with xt/d-°-f- = 0.64/62. The streak 
exposure time correction is found to be ~ 240, suggesting 
that the true count rate of M82 X-l is 2.1 c/s, and its 
un-absorbed luminosity is 8.7 x 10 40 erg s -1 . 

The spectrum of SI has NH$i = 1.081q 24 x 10 22 cm -2 
and Tsi = 1 . 12^21 with xt/d-o.f. = 1.38/48. Placed at 
the distance of M82, SI would have an un-absorbed lu¬ 
minosity of 7.2 x 10 38 erg s _1 , ~ 1% of that of M82 X-l. 
However, due to the severe pile-up of M82 X-l, the 0.3-8 
keV count rate of SI is almost half of M82 X-l’s count 
rate, enough to cause a significant positional displacement 
of M82 X-l. 

To test the spectral difference for subregions of the PSF, 
we have extracted the spectra from two subregions of the 
nearby bright source J095551.2 in parallel to X-l and SI 
of M82 X-l (see Fig[l]). As expected, these spectra are 
similar to each other with T ~ 2.0. 

As shown above, Si’s spectrum is consistent with those 
in the low-hard states of black hole X-ray binaries (Remil- 
lard & McClintock 2006) and is intrinsically different from 
that of M82 X-l or its readout streak at the 90% confidence 
level, again supporting that SI is a real source. 

2.3. Temporal Behavior of SI 

We further analyze other Chandra ACIS observations in 
a similar way as in Obs ID 10543. SI is not detected in any 
other ACIS observations, including the three observations 
10542, 10925 and 10544, observed about a week before and 
after Obs ID 10543, respectively. Using the same 0.7" ex¬ 
traction region and the same annular background region, 
we extract counts at Si’s position for all other ACIS ob¬ 
servations. In Obs IDs 378, 379, 380, 6097, 8190, 10025, 
10026, 10027 and 10545, SI was fully buried in the PSF 
wings of X-l, so extraction results from these observations 
are excluded. In other observations, the PSF wing of M82 
X-l still partly cover SI regions, and these obvious X-l 
PSF wing regions are excluded. After this process, it is 
still possible that some photons in the extraction regions 
could come from X-l, thus the estimated count rates of SI 
should be treated as upper limits. In Fig[2]we present the 
results. The 0.3-8 keV count rate of SI was < 0.001 c/s in 
1999 (Obs ID 361) and was < 0.007 c/s in 2009 (Obs ID 
10542), then rose to ~ 0.016 c/s within a week (Obs ID 
10543) and dropped to <0.0007 c/s with 7 days (Obs ID 
10544). Judged from the light curve, SI is likely a tran¬ 
sient source in 0.3-8 keV band. This transient behavior, as 
well as the fact that SI was usually in the shadow of M82 
X-l’s PSF is likely the reason why SI was only detected 
once in 22 Chandra observations. 

2.4. Astrometry 


4 see http://cxc.harvard.edu/cal/Hrc/PSF/acis_psf_2010oct.html and http://cxc.cfa.harvard.edu/ciao/caveats/psf_artifact.html for details 
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It is now possible to identify the radio flare source given 
the centroids of M82 X-l and SI. Using three bright ACIS 
point sources (a > 60; see Table [2] for details) with op¬ 
tical jnahchesfrom the Sloan Digital Sky Survey(SDSS, 

I Aim et aJ1l2014f) . we find a ~ 0.4" displacement between 
the ACIS coordinate frame and the SDSS coordinate frame 
(the same as the VLBI coordinate frame with an error of 
0.2"). The position of the radio flare source is corrected ac¬ 
cordingly and registered to the ACIS image (Fig [I]). The 
flare source is now 0.96" away from M82 X-l, but only 
0.28" away from SI. Clearly, the position of M82 X-l is 
no longer consistent with that of the radio flare source. 
On the other hand, Si’s position is consistent with the ra¬ 
dio flare source within the positional uncertainty (~ 0.4"; 
mainly from the astrometric correction). 

If the radio flare is indeed associated with SI, it can be 
the radio outburst from a black hole X-ray bi nary w hen the 
binar y transits f rom hard to soft states (e.g., iFender et al.l 
l2004lh Recently. iMiddleton et al.l (1 201 3li discovered strong 
radio emission at a flux level < lrnJy from a micro-quasar 
in M31 during its X-ray outburst at an X-ray luminosity 
of 1.3 x 10 39 erg s _1 . The properties of SI are consistent 
with those of the M31 micro-quasar, and it can well be a 
similar micro-quasar, albeit with higher radio flux (9mJy). 
However, unlike in the case of the M31 micro-quasar, the 
radio flare and the X-ray outburst for SI were not detected 
simultaneously, and further simultaneous monitoring ob¬ 
servations in radio and X-ray are needed to determine its 
exact nature. 

3. DISCUSSION 

M82 X-l, one of the most promising IMBH candidates 
in the local universe, is possibly associated with a radio 
flare source, suggestive of relativistic beaming effects for 
its extremely high X-ray luminosities. We have applied 
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the sub-pixeling technique to the longest Chandra obser¬ 
vation of M82 X-l, and split it into two separate discrete 
sources, namely X-l itself and SI. Both relative positions 
and spectral properties show that SI is a real source and 
not a PSF artifact, although it was detected only once 
out of 22 Chandra observations. Careful astrometric work 
shows that the radio flare source is associated with the 
transient SI. M82 X-l is no longer associated with the ra¬ 
dio flare, eliminating the necessity for relativistic beaming. 
The radio flare and X-ray outburst of SI are analogous to 
those of a recently discovered micro-quasar in M31 (Mid¬ 
dleton et al. 2013), making it a likely micro-quasar in the 
shadow of M82 X-l. 

The discovery of SI, however, does not affect previous 
results on M82 X-l’s spectral properties significantly (e.g., 
Feng et al. 2010). This is because SI luminosity, even 
in its outburst, is only ~ 7 x 10 38 erg s _1 . In compari¬ 
son, M82 X-l exhibited X-ray luminosities from 2 x 10 4 ° 
to > 10 41 erg s _1 (e.g., Kaaret et al. 2006; Feng et al. 
2010), 30-100 times higher than the outburst luminosity 
of SI. Furthermore, SI was detected only once out of 22 
Chandra observations between 1999 and 2014, and it must 
remain in a quiescent state at even lower luminosities most 
(>95%) of the time . However, SI could become impor¬ 
tant for on-axis observations when M82 X-l suffers severe 
pile-up. In the case of Obs ID 10543, the count rate of SI 
is almost half of the observed count rate of M82 X-l. 
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Fig. 1.— 1/8 pixel Chandra 0.3-8 keV image of M82 X-l. Solid black circles are 3-cr error circles for detected sources as noted. Solid green 
circles, dashed green annulus and rectangle are for source and background spectra extraction. The Chandra instruments +Z and readout 
streak directions are given at the upper-left corner. The cyan region shows the possible spurious source region caused by the PSF distortion. 
The ‘X’ and the dashed black circle mark the radio flare position and its uncertainty brought by astrometry correction, respectively. The two 
sources J095551.2 and J095550.6 marked with solid black circles at the lower-left corner are two nearby bright X-ray sources. Note that there 
are no count excess in the lower left direction of J095551.2. 



5.1X10 4 5.2x10 4 5.5x10 4 5.4x10 4 5.5x10 4 5.6x10 4 5.7x10 4 

Observation Date(MJD) 

Fig. 2.— Averaged 0.3-8 keV net count rates of SI in Chandra ACIS observations between September 1999 to February 2013. The 12 
points from left to right refer to Obs ID 361, 1302, 2933, 5644, 6361, 10542, 10543, 10925, 10544, 11104, 13796 and 15616, respectively. The 
filled triangle corresponds to the detection of SI in Obs ID 10543, and the opened triangles are from observations without detection and 
should be treated as upper limits. 
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Table 1 

Results of spatial and spectral analysis of M82 ULX-1, the readout streak and SI in 0.3-8 keV. An absorbed 

POWER LAW MODEL IS USED TO CHARACTERIZE THE SPECTRA. 


Source 

RA 

Dec 

Net Counts 

Significance 

NH 

r 

X 2 /d.o.f. 






(1.0 x 10 22 cm -2 ) 



ULX-1 

9:55:50.123 

+69:40:46.54 

4165 

238a 

0.56 

-0.68 

2.1/174 

ULX-1 Streak 

- 

- 

- 

- 

1 oe:+0-38 
i.zo_o.33 

1 8 C ) +0,35 
±.OO_ 0 33 

0.64/62 

SI 

9:55:50.245 

+69:40:45.67 

2190 

71a 


i i O+0.22 
1 -- l ^-0.21 

1.38/48 


Table 2 

Positions of sources for ACIS-HRC and ACIS-SDSS astrometric registration. 



Selected Sources Detected in Both ACIS and HRC Observations 


ACIS RA 

ACIS Dec Net Counts Significance HRC RA HRC Dec 

Displacement 

-r, nail - 


9:55:50.689 
9:55:51.257 
9:55:50.355 
9:55:51.479 
9:55:47.488 
9:55:14.583 b 


+ 69:40:43.61 
+ 69:40:43.79 
+ 69:40:36.36 
+ 69:40:35.90 
+ 69:40:59.61 
+ 69:47:35.55 


1744 

7063 

1747 

2252 

587 

2201 


103cr 
501 c 
86c 
105c 
61er 
131c 


9:55:50.803 

9:55:51.395 

9:55:50.473 

9:55:51.591 

9:55:47.585 

9:55:14.496 


+69:40:43.37 

+69:40:43.58 

+69:40:36.18 

+69:40:35.70 

+69:40:59.42 

+69:47:35.34 


0.63 

0.73" 

0.63" 

0.60" 

0.53" 

0.49" 



Bright ACIS X-ray Sources identified 

in SDSS Observations 


ACIS RA 

ACIS Dec 

Net Counts 

Significance 

SDSS RA 

SDSS Dec 

Displacement 

9:56:58.658 

+ 69:38:52.13 

652 

90cr 

9:56:58.678 

+69:38:52.54 

0.422" 

9:55:05.206 

+69:44:42.47 

574 

61c 

9:55:05.227 

+69:44:42.85 

0.394" 

9:55:14.583 b 

+69:47:35.55 

2201 

131c 

9:55:14.546 

+69:47:35.85 

0.354" 


Note. — a M82 ULX-1. 

b The source detected in HRC, ACIS and SDSS observations. 















